Tracking Quasiparticle Energies in Graphene with Near Field Optics 
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Advances in infrared nanoscopy have enabled access to the finite momentum optical conductivity 
cr(q, uj). The finite momentum optical conductivity in graphene has a peak at the Dirac fermion 
quasiparticle energy e(&F — q), i.e. at the Fermi momentum minus the incident photon momentum. 
We find that the peak remains robust even at finite temperature as well as with residual scattering. 
It can be used to trace out the fermion dispersion curves. However, this effect depends strongly on 
the linearity of the Dirac dispersion. Should the Dirac fermions acquire a mass, the peak in a(q, w) 
shifts to lower energies and broadens as optical spectral weight is redistributed over an energy 
range of the order of the mass gap energy. Even in this case structures remain in the conductivity 
which can be used to describe the excitation spectrum. By contrast, in graphene strained along 
the armchair direction, the peak remains intact, but shifts to a lower value of q determined by the 
anisotropy induced by the deformation. 



I. INTRODUCTION 

Graphene, first isolated in 2004^ has been the host of 
a variety of novel electronic properties. The main dif- 
ference in graphene is its unique energy dispersion. The 
charge carriers in graphene are massless Dirac fermions, 
which accounts for the differences from the conventional 
2D electron gas. Remarkable behaviour has already been 
reported in the plasmon dispersion relation^ - — as well as 
the optical conductivity,— ~~ which supports a transverse 
electromagnetic mode3 There is also recent evidence for 
plasmarons ; 10 ' 1 1 a new type of quasiparticle formed by 
the interaction of charge carriers with plasmons. Opti- 
cal spectroscopy is a useful tool for obtaining information 
about the dynamics of charge carriers, and has been used 
to great success in graphene^ 

The real part of the q —¥ optical conductivity in 
graphene is well known. At finite chemical potential, \x, 
it contains a Drude peak at uj = due to intraband ab- 
sorption, followed by a Pauli-blocked region. There is 
then a sharp rise at uj = 2/i to a universal background 
conductivity (To = e 2 /4^^— due to interband transi- 
tions. Experimentally, the region which should be Pauli- 
blocked and have no absorption does not fall below about 
(To /3 . 1(> Electron-electron interactions , electron-phonon 
interactions, and impurity scattering can all provide con- 
tributions to the optical conductivity in this region, but 
nothing as large as the observed value J2r— The origin of 
this anomalously large background is still a mystery. 

More recent experiments have granted access to 
the finite momentum transfer optical conductivity 
cr(q, uj)^t— In the paper by Fei et ali^i they describe 
how an atomic force microscope (AFM) operating in tap- 
ping mode allows one to obtain information about the fi- 
nite q conductivity. The incident light scatters off the tip 
and is confined to a nanoscale region. The precise details 
depend on the geometry of the tip, and Fei et al. report 



a distribution of q values with a peak at q ~ 3.4 x 10 s 
cm . In principle, a sharper tip would lead to higher 
confinement, and thus larger values of q could be acces- 
sible through adjustments to the AFM tip. 

In this paper we study the properties of the quasipar- 
ticle peak at uj = q in the real part of the optical con- 
ductivity. In section [II] we discuss the properties of this 
peak and its relationship to the joint density of states at 
T = as well as at finite temperatures. We find that the 
quasiparticle peak remains robust even at large tempera- 
tures. In section HIH we consider the presence of residual 
scattering and provide a simple analytical formula for the 
quasiparticle peak. Again, the peak position remains ro- 
bust, even for large impurity concentration. We finally 
consider two methods of altering the Dirac dispersion in 
section HVl gapped and strained graphene. We find that 
these modifications do alter the position of the quasipar- 
ticle peak and so, the linearity and isotropy are crucial 
for its robustness. 



II. FORMALISM AND EXPRESSIONS FOR 
THE CONDUCTIVITY 

The xx component of the real part of the finite tem- 
perature optical conductivity is given by 



cr xx {q,uj) 



lf(uj' + uj)-f(uj')}(hj' 



~2T 



F ss > (4)A S (k, uj')A s ' (k + q, uj' 



(1) 

In the above, Co = e 2 /4 is the universal background con- 
ductivity, F ss i ((f>) are the coherence (or chirality) factors, 
f(uj) = l/(e l3u> + I) is the Fermi-Dirac distribution func- 
tion, and A s (k, uj) is the spectral density. We chose our 
x-direction along the zig-zag axis (see Figured]). From 
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FIG. 1: (Color online) (top) The scattering geometry in k- 
space for the optical conductivity, and polarization. Both 
the optical conductivity and polarization contain coherence 
factors F ss i((f>) = i(l + ss' cos <j>). The important difference is 
that <j) — 6k + 6k+q in the optical conductivity, while <j) — 6^ — 
Ok+q in the polarization. We orient our axes so that x is along 
the so-called zigzag direction of graphene, while the i/-axis is 
along the armchair direction, (bottom) Energy dispersions 
for bare, strained and gapped graphene. The effect of strain 
distorts the Dirac cones from a circular to an elliptical cross 
section. Gapped graphene retains its circular shape, but the 
low energy dispersion is now quadratic, and the Dirac point 
is no longer accessible. 



here on we use <r(q, uj) to denote a xx (q,uj). We work in 
units where h = vp = 1. 

The spectral densities, A s (k,w), reduce to Dirac-delta 
functions in the bare band case. In the presence of inter- 
actions described by a self energy E s (k, ui) they are given 
by 



where et 
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|Im£ s (k,w) 



7T (u> - Re£ s (k,w) - eg.) 2 + |Im£ 8 (k,w)| 



(2) 



The peaks in the spectral density control the shape of 
the optical conductivity and correspond broadly to two 
types of processes: intraband and interband transitions. 
The intraband transitions occur at uj = q and are the fo- 
cus of this paper. Interband scattering is responsible for 
subsequent peaks in the spectral functions which occur 
at ui = 2/i — q and uj = 2/i + q. We will first consider 
the bare band case, and examine the effect of impurity 
scattering in Section Mil 



A. Results for Bare Bands, T=0 

Since the spectral functions are simply given by Dirac- 
delta functions in the bare band case, the physics is gov- 
erned by the coherence factor F ss > (</>), which encodes in- 
formation about scattering. It is given by 



F ss >(4>) = i(l + ss' cos < 



(3) 



The angle <f> is defined in terms of the angles of k and 
k + q, denoted by #k and 6>k+ q respectively (Figure [T|) . 
For the optical conductivity <p — $k + #k+q, in contrast 
to the polarization, in which <f> = #k — #k+q- 

We can write down an expression for F ss > (</>) in terms of 
the magnitudes of k and q and their angles with respect 
to the fc^-axis, 9 and a respectively (Figured]). We have 



1 / , fccos(2g)+gcos(fl + a) \ 

^ 1 + ss — ; i (4) 

2\ ^k 2 +q 2 + 2kqcos{9 -a) I 



for the optical conductivity, and 



F.A<t>) 



k + qcos(9 — a) 



^/k 2 + q 2 + 2kq cos(6> - a) 



(5) 



for the polarization, which we include for comparison. 
Notice that the polarization only involves the angle 9 — a 
while the optics contains 9 — a, 9 + a, and 29. A con- 
sequence of this is that for the polarization, the depen- 
dence on the direction of the photon momentum q can 
be eliminated by a shift in the integration variable. For 
the optical conductivity no simple change of variables ex- 
ists, and F ss i ((/)) remains dependent on the angle q makes 
with respect to the k x -&xis. 

In an isotropic system only two directions need to be 
considered for q: q along k x (longitudinal) and q along 
k y (transverse). This gives 



F 8 A4>) = \[i 



, kcos{29) + qcos(9) 

TFT 



2kqcos(9) 



for the longitudinal part of a xx , and 



F ss .{$) = - 1 



, kcos{29) - qsm(9) 
^k 2 + q 2 + 2kqsin(9) 



(6) 



(7) 



for its transverse part. The difference between Eq. [6]and 
Eq. [7] has a drastic difference in the shape of cr(q, u>), its 
longitudinal part diverges at uj — q, while the transverse 
parts vanishes at ui = q (Figure [5]) • 

To understand this difference in behaviour we intro- 
duce a reference function, the joint density of states 
(JDOS). It is given by 



JDOS(q,cj) = 4 



d 2 k 



- m) - JK+q - m)] 

x5(w + e k -e k ' +q ), (8) 
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FIG. 2: (Color online) The real part of the longitudinal and 
transverse conductivity, a L and a T , for q = 0.4 scaled by 
q/ir as a function of cj//t. JDOS and the Polarization II as a 
function of w//t. The pre- factors have been chosen to make 
the longitudinal conductivity and JDOS agree at cj = q. 



where the factor of 4 is from the degeneracy (spin and 
valley) in graphene. The JDOS for the intraband piece 
at T = is given by 



JDOS(q,u;) = 



1 



Q(w -q+ 2/i) f (2/i + w)y/(2fi + w) 



2_ q 2 



_ {q2 2w2) ln yszM 



the final state to determine the contribution to the con- 
ductivity. For the longitudinal case, that is, q taken along 
the k x direction, the final state momentum is the sum of 
the magnitude of the initial momentum k and the pho- 
ton momentum q and results in a maximum momentum 
along k x . For the transverse case, the initial k must also 
be transverse to be an allowed transition (recall u> = q). 
The resulting state has no-momentum along the k x di- 
rection, and so the transverse conductivity vanishes. 

The joint density of states, Eq. [SJ bares a strong 
resemblance to both the polarization II, and the 
longitudinal conductivity, which have been computed 
previously^^^ The conductivity, polarization and 
joint density of states are all shown in Fig. [5J The factor 
q/ir multiplying the conductivity was chosen to make the 
functions agree at oj = q, and the factor multiply- 
ing the polarization was chosen to make the polarization 
have the same prefactor as the JDOS. 

The agreement between the joint density of states and 
the polarization is excellent (Figure [5]). In fact, the dif- 
ference is given by 



JDOS - -n =- 

71" 2^ 2 ^Jlf 



hi- 



IV 



+ 2fi+ y/(2n + w) 2 -q 2 



2/i - w + \J(2\x - w) 2 - q 2 



(10) 



and so the differences between the two are logarithmi- 
cally small. In Fig. 3 we also see that ^JDOS agrees 

remarkably well with ^ a L . The difference between these 
two functions is in fact controlled by the same logarith- 
mic factor as in Eq. [TUJ 



- 8(2/1 — q — w) f (2/i - w) v / (2/i-w) 2 - q 2 
^/(2[l-w) 2 -q 2 + (2/i - to) N 



-(<T - 2«T)ln 



q 



(9) 



The longitudinal and transverse conductivity differ 
from the joint density of states by the coherence factors, 
as mentioned above. Both the joint density of states and 
the conductivity contain the same delta function. Eval- 
uating the coherence factors Eq. [S] and Eq. [7] at lu = q 
subject to the delta function constraint shows that F ++ 
= 1 and respectively. This explains the fact that the 
longitudinal conductivity has a square root singularity 
(inherited from the JDOS), while the transverse conduc- 
tivity is zero. 

A physical picture for this difference is as follows. We 
consider a possible optical transition with momentum 
transfer q and energy u = q from an occupied state be- 
low the chemical potential, to an empty state above the 
chemical potential. We will consider the momentum of 



B. Finite Temperature 

We now turn to the effect of finite temperature. In this 
case, the expression for the real part of the longitudinal 
conductivity is 



o- L (q,w) 



(To 



— T / d 2 k[/(4- M )-/(4; 

7TLJ £ — 4 J 

ss' 

x F ss ,{(t))5{Lo + e{ 



-m)] 

s k+q); 
(11) 



where F ss i{4>) is given as in Eq. [6] We use the delta- 
function to do the integral over the angular variables, 
and find that the conductivity naturally separates into 2 
parts, one part for u < q and the other for u > q. They 
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are given by 

8 w 
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dx 



sinh ^ 



cosh + cosh 2±|^ 



rZ.r 



sinh 



cosh $ + cosh 9+2 f+ 2 ^ 



\/x(a; + g) 



if 



7T q 2 ^Jq 2 - w 2 



dx 



sinh 



2T 



COSh £ + COSh 2± 2 ^2M 



y/x(x + (?) 

(12) 



for uj < q, and 
o> 8 



^ q 2 \Jv? 



dx 



sinh 



2T 



cosh $ + cosh ? - 2 2 3 r - 2A ' 



(13) 

for w > q. We have simplified the expression for ct< by 
noting that the thermal factors in the second term cause 
it to be much smaller than the first. 

The remaining integrals were evaluated numerically 
and the results for T/fx = 0, 0.03, 0.07, 0.15, and 0.3 
are shown for momentum transfer q/kp = 0.4 in Fig. |3] 
There is a sharp quasiparticle peak from intraband tran- 
sitions at u) = q. It remains sharp even at elevated tem- 
peratures. However, the interband transitions are ther- 
mally broadened. One naively expects thermal broad- 
ening to occur over a width ~ T, and we see the effect 
is much larger for the interband transitions. This ex- 
cess broadening can be understood as an enhancement 
from the square root singularity present in the JDOS. 
All the finite temperature curves intersect at the point 
lu = 2/i. This could be used as a method of determining 
the chemical potential. We also computed the optical 
spectral weight (see insets Fig. |3J) given by 



I(oj) = doJ 
Jo 



,cr(q,w') 



CO 



(14) 



We see that finite temperature shifts spectral weight from 
the interband region into the previously forbidden region 
q < ui < 2/i — uj. Fig. [3] also shows the finite tempera- 
ture effect for momentum transfer q/kt = 1.0. Again the 
quasiparticle peak is sharp, and the interband transitions 
are smeared. Notice that for larger q the spectral weight 
carried by the quasiparticle peak is diminished. Increas- 
ing q has the effect of decreasing the spectral weight car- 
ried by the quasiparticle peak. This spectral weight is 
regained in the interband transitions so that the optical 
sum rule remains satisfied. 




FIG. 3: (Colour online) (Top panel): The real part of the 
finite momentum optical conductivity <r(q, ui) as a function 
of u/n for q/k F = 0.4 and T/fi = 0, 0.03, 0.07, 0.15, 0.3 
for bare bands. There is a strong quasi-particle peak at u = 
q, which is unaffected by the finite temperature. The finite 
temperature smears the interband contribution, and begins to 
fill in the Pauli-blocked region for large enough temperature. 
(Bottom panel) : The real part of the finite momentum optical 
conductivity a(q, w) as a function of uj/fi for q/kF = 1.0 and 
T/fi = 0, 0.03, 0.3 for bare bands. Included for comparison is 
the T = result with a residual scattering rate = 0.005. 
(Insets): The insets show the optical spectral weight for T/fi 
= 0, 0.3. The quasi-particle peak holds less spectral weight at 
larger momentum transfer (q). For the q/kF = 0.4 case, we 
can see that the finite temperature has transferred spectral 
weight to the previously Pauli-blocked region. 



5 



Although we cannot obtain an analytic formula for the 
spectral weight carried by the peak at general q, we can 
obtain expressions for the spectral weight at T — in 
the limit q — > 0, for both the longitudinal and transverse 
conductivity. For the quasiparticle peak, in the q — > 
limit we obtain 



9 z „ =?> 

Co nq^y'q^ — w z 

— (16) 

So that 

/ cfo; — = / dw — = 2fi, (17) 
Jo a o Jo c 

and both the transverse and longitudinal peaks carry the 
same spectral weight in this limit. The fact that the 
quasiparticle peak in the longitudinal and transverse con- 
ductivities have the same spectral weight, combined with 
isotropy implies that the peak carries with weight regard- 
less of the direction of q in the q — > limit. 



III. EFFECT OF IMPURITIES 




FIG. 4: (Color online) The real part of the longitudinal optical 
conductivity a(q,(jj) for q/kf = 0.4 as a function of cj/fi. 
Included are the result for bare bands at T = in black, a 
numerical evaluation of Eq. [T] including impurity scattering 
with j/fi = 0.005 in dashed green, and our analytic expression 
for the quasiparticle peak, Eq. [5S]in light green circles. 



We have, until now, considered only the bare band 
case. We now consider the effect of scattering. The sim- 
plest approximation which includes scattering is to take a 
self energy with ImS(k, oj) — 7 and Re£(k,u;) = 0. Note 
that, in particular, we ignore vertex corrections. In this 
case, and at T — 0, the general formula for the intraband 
conductivity, Eq. [TJ becomes 



— = — f doj' [ kdk [ — 
Co w J_ u J Jo 2tt 



x F ++ (0) 



1 



7 



7T 2 7 2 + (u' + k + fi) 2 

7 



7 2 + (u> + uj + y/J? + q 2 + 2kqcos9 + ( u) 2 

(18) 

We are only interested in the case when uj,j,q <C fi. 
In this case, consideration of the Lorenzian factors, tells 
us that the dominant part of the integral in Eq. [T5] is 
from the region k « fi. This allows us to simplify the 
expressions for both F++ as well as the second Lorentzian 
in Eq. [18] Working to lowest order in q we obtain 

Co w J-uj J Jo 2tt 



x cos 2 (6>) 



1 



7 



7T 2 7 2 + (w> + k + /i) 2 

7 



7 2 + (u)' + uj + k + /i) 2 ' 



(19) 



where uj = uj + qcos{9). The integration over k can be 
performed and we have 
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dw' 



cr n z w J_ w J 2n 



d6 



1 



w 2^j 2 + w 2 
{uj' + fi) 2 + 7 2 



{uj' + uj + fi) 2 + 7 2 

6 / / - \ _i ( w' + Q 

47 {uj + uj + fi) tan h fi 

7 V 7 



+ -(u/ + /z)tan 1 ( — + M 
7 V 7 



(20) 



Under our conditions that w,7, q <C fi this simplifies 



to 



a L 4/i 



<m- 



2 7 cos 2 6) 



Co I -2 Jo {w + q cos 8) 2 + 47 2 
A similar calculation gives the transverse conductivity 

—T a,, /-27T 1„,„«~2/1 



(21) 



- = ^ / dfl- 



27 cos 



°"o 7T 2 Jo { w + <zsin#) 2 + 47 2 ' 
and the Polarization 

f27T 



n 



UJfL 



de 



20 



7T 2 7 ( w + qcosB) 2 + Aj 2 



(22) 



(23) 



These integrals can all easily be evaluated. They are 
most conveniently expressed in terms of the complex 
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FIG. 5: (Color online) The quasiparticle peak in the real part 
of the optical conductivity for momentum transfer q = 0.4 
(main panel) and q — 1.0 (inset) as a function of u/fi. We 
show impurity concentrations of j/fj, = 0, 0.05, and 0.005. 
The position of the peak stays robust even with large impurity 
content and is broadened as the impurity content increases. 



number 



y/(w — q + 2i r y)(w + q + 2ij) 



(24) 



Our final expressions for the quasiparticle peak in the 
presence of scattering are 



a 



I67/1 



1 - 2wRc(Z) + —Im(Z*) + 2- f lm(Z) 
2 7 



a 

CO 



I6/X7 

irq z 



Jm(Z*) 
2tW 



1 



l m (n) = ^Im(Z*). 

IT 



(25) 
(26) 

(27) 



The agreement between the analytic expressions and a 
numerical calculation of Eq. [T]is excellent (Fig. 0|). We 
verified that this agreement is maintained up to q = 1.0, 
even though our derivation assumed q was a small param- 
eter. We show evaluations of the longitudinal conductiv- 
ity, Eq. 1251 for two impurity concentrations in Fig. [5] for 
q =0.4 and 1.0. The case with no impurity scattering is 
included for reference. The peak becomes progressively 
broadened as we increase the impurity scattering rate, 
7, but the peak position remains robust, even for large 
disorder. 

Interestingly, our formulas for the conductivity with 
finite residual scattering are almost the same as the q — > 
limit of the conductivity given in Eq. [15] and Ql)] with 
the replacement u — > uj + 2*7. There is an additional 



term present in our formulas that is not captured by this 
simple substitution. For the longitudinal conductivity 
our expression is the same as (compare with I15[) 



^ = 8 M Re( ^ + 2ll? 



(28) 



While for the transverse conductivity it is the same as 
(compare with [T| 



— = 8^Re 
Co 



q 2 — (bj + 2i~f) 2 \ I67/Z 



nq z 



irq z 



(29) 



As a final remark, we comment on the difference be- 
tween the polarization and the optical conductivity. In 
the non interacting case it has been shown^ that the 
polarization is related to the conductivity through the 
standard formula 



^jim n . 



(30) 



We remarked on the differences between the coherence 
factors of the polarization and the conductivity in section 
III Al In the non-interacting case, the spectral densities in 
Eq. [1] reduce to delta functions, and the delta function 
constraints restrict the coherence factors of the polariza- 
tion and the conductivity to be related through only the 
factor ui/q 2 . In the presence of impurities, the delta func- 
tions become broadened and the coherence factors are no 
longer proportional. In fact, we see that the replacement 
given by Eq. [3D] using the polarization in the presence of 
impurities Eq. [27]only generates one of the terms present 
in the optical conductivity Eq. 1251 It is worth noting that 
the term generated by the polarizability is the dominant 
term near u = q. At small values of u the Z independent 
piece becomes the dominant contribution. All the terms 
not proportional to the polarizability are suppressed by 
factors of 7 so that the correct limit is obtained as we 
turn off impurity scattering. 



IV. MODIFICATIONS TO THE DIRAC 
SPECTRUM 

Finally, we examine the consequences of altering the 
energy spectrum in graphene on the quasiparticle peak 
in the optical conductivity. We will consider two phys- 
ical mechanisms for altering the spectrum in graphene. 
The first is the opening of a mass gap, A. The second 
is the application of strain, which makes the Fermi ve- 
locities along x and y different. These alterations to the 
spectrum are shown pictorially in Fig. [T] 



A. Gapped Graphene 

In graphene with gap A, the energy eigenvalues are no 
longer linear in k but instead are given by 



= y/k 2 + A 2 



(31) 
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FIG. 6: (Color online) The quasiparticle peak in the real part 
of the optical conductivity as a function of uj/h for q/kF = 0.4 
for graphene with mass gap A = 0, 0.05, 0.1, 0.15, 0.2, and 
0.25. The peak is pulled to smaller values of q and slowly 
broadened as the mass gap increases. The broadening hap- 
pens over an energy scale approximately given by the mass 
gap, and is physically caused by changes in the joint density 
of states (shown in the inset). 



The optical conductivity of gapped graphene was first 
studied by Scholz and Schliemamii^i, As we saw in section 
HT1 the peak in the optical conductivity depended strongly 
on the joint density of states. The joint density of states 
for gapped graphene is 



JDOS(q,cj) = 



1 



4-7T 2 \J q 2 — w 2 
Q(w - qx + 2/i) \ {2fi + w)^{2n + w) 2 - q 2 x 2 (32) 



(q 2 xl - 2w 2 )ln 



^/(2^ + w) 2 -q 2 x\ + (2/i + w) 



qx 



-6(2/x — qx — w) M2/i - w)yj (2/i - w) 2 — q 2 x\ 
-tfxl - 2w 2 ) In V(2^-H 2 -^ + (2M- W ) ^ 



where /i = \Jk 2 F + A 2 and xq = J\ 



4A 2 



qx Q 

'(33) 
Figure [5] 

shows the quasiparticle peak in the optical conductivity 
for several values of the gap, as well as the joint density 
of states. We see that as the gap opens, the joint den- 
sity of states flattens out and is pulled back to smaller 
values of uj. Consequently, this behaviour is inherited in 
the optical conductivity. The peak is shifted to smaller 
values of uj and broadened as A increases. In particu- 



lar, the flattening onsets at uji = 
and persists to uj u = e{kp + q) 
conductivity then vanishes. 



A 2 



- y/k% + A 2 



- e(k F - q) 
where the 



B. Strained Graphene 

We consider, for simplicity, the case where strain is 
applied along the armchair (or y) direction in graphene 
(Fig. [1]). The effect of such a strain can be captured 
by introducing two strain parameters, ^ x and "f y , which 
control the anisotropy into the Dirac Hamiltoniart 2 ^— 



^yOyky. 



(34) 



First we consider the longitudinal conductivity, that 
is, q along k x . By examining the Kubo formula for the 
conductivity, Eq. [TJ we find that the change of vari- 
ables k = {^xkx,lyk y ) is sufficient to give us a result for 
the strained conductivity. There is a Jacobian from the 
fc-space integration which contributes a factor 1/(7^7^). 
Additionally, there is a factor of v 2 F (1 in our units), which 
contributes a j 2 (for u xx ). Lastly, since q appears only 
in 6k+ q it is changed to q = q"/ x . Expressed in terms of 
k,q all the integrations are the same as the unstrained 
case. Thus, we arrive at the simple result 



strained 



— 0teo(?7s,w). 

7y 



(35) 



A similar calculation for the transverse conductivity gives 



strained 



few) 



lx 
— ( 



(36) 



We see that the position of the peak is shifted from oj = q 
to uj = j x q, and the overall conductivity is modified by 
the ratio ^ x /ly 

This shift in the peak is also understandable from the 
physical picture described in section III Al described ear- 
lier. Considering the same transitions described there 
gives the peak at uj — q. The effect of strain on the sys- 
tem is to distort the shape of the cone. Focusing on the 
x direction, we have that, geometrically, this changes the 
lengths of the vectors by a factor j x . To still land on the 
energy dispersion (and thus be an allowed transition), 
the energy must also be modified by a factor of j x . This 
simple geometric consideration gives the shift in the peak 
position. 



V. CONCLUSIONS 

We have considered the peak in the real part of the 
near-field optical conductivity. The peak is located at 
uj = q, and, as long as the dispersion is linear and 
isotropic this position is robust. This quasiparticle peak 
is due to intraband transitions, and is the finite q ana- 
logue of the Drude peak, in the q — conductivity. At 
q = the peak carries 2/i worth of spectral weight so 
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that the optical sum rule is satisfied. At finite q the 
quasipatricle peak carries less weight, with the missing 
weight transferred to interband processes. 

We find that both finite temperatures and finite resid- 
ual scattering rate fill in the Pauli-blocked region, and 
that this filling is enhanced near u) = q due to a square 
root singularity in the density of states. We used the 
Kubo formula in the bubble approximation, which ig- 
nores vertex corrections, to understand how the bare 
band picture is modified, in the presence of residual scat- 
tering. This allowed us to derive simple expressions for 
quasiparticle peaks in longitudinal and transverse con- 
ductivity, as well as the contribution from intraband scat- 
tering to the polarization. 

Finally, we examined the effect of altering the Dirac 
dispersion on this peak. We found that in the presence 



of a gap, the peak was shifted to smaller values of u, 
and reduced in size, no longer able to feel the square 
root singularity in the JDOS. The peak also broadened 
as the mass gap increased. In the presence of strain, the 
effect was to scale the height of the peak, by a geometric 
factor related to the ratio of the anisotropics induced by 
the strain, in addition to the position of the peak being 
shifted by the relevant anisotropy parameter. 
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